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ABSTRACT 

We present cosmological hydrodynamical simulations of the formation of dwarf galax- 
ies in a representative sample of haloes extracted from the Millcnnium-II Simulation. 
Our six haloes have a z = mass of ~ 1O 1O M0 and show different mass assembly 
histories which are reflected in different star formation histories. We find final stellar 
masses in the range 5 x 10 7 — 10 8 Mq, consistent with other published simulations of 
galaxy formation in similar mass haloes. Our final objects have structures and stellar 
populations consistent with observed dwarf galaxies. However, in a ACDM universe, 
1O 1O M0 haloes must typically contain galaxies with much lower stellar mass than our 
simulated objects if they are to match observed galaxy abundances. The dwarf galaxies 
formed in our own and all other current hydrodynamical simulations are more than an 
order of magnitude more luminous than expected for haloes of this mass. We discuss 
the significance and possible implications of this result. 

Key words: cosmology: theory - galaxies: dwarf - galaxies: formation - galaxies: 
evolution - galaxies: luminosity function, mass function - methods: N-body simula- 
tions 



1 INTRODUCTION 

Dwarf galaxies are by far the most abundant type of galaxy 
in the Local Group and in the Universe. They span a large 
range of stellar masses, morphologies and star formation his- 
tories. The largest dwarf irregulars such as the large Mag- 
ellanic Cloud have stellar masses of ~ 10 9 Mo, rotationally 
supported and Hi-rich disks, and strong ongoing star for- 
mation. In contrast, dwarf spheroidal galaxies have stellar 
masses from 1O 7 M0 to below 1O 3 M0, they possess no inter- 
stellar gas, and they show no sign of rotational support or 
ongoing star formation. 

The number of dwarf galaxies observed in the Local 
Group continues to gro w as new, 'ultra-faint' satellite galax- 
ies a r e discovered (e.g . iMartin et al.l 120061 : IChapman et al.l 
120071 : iBelokurov et all |20ld h Estimates using luminosity 
functions corrected for completeness and bias predict the 
total numb er of faint satellites to be an order of mag nitude 
higher still (|Tollerud et al.ll2008l ; IKopqsqv et al.ll2008r ). Nev- 
ertheless, this is still much smaller than the total number 
of dark matter subhaloes found in high-resolution Simula- 



tions of the standard ACDM cosmoloe 


y (e.g. iKlvpin et al. 


199S 


:lMoore et al.lll999; Dicmand et al. 


2007;ISprhi£el et al. 


2008 


). This difference has become known as the "Missing 



Satellites Problem". It may only be an apparent discrep- 
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ancy, however, since it can be removed if one accounts for 
the fact that not all low-mass subhaloes must contain stars, 
and those that do may have very high mass-to-light ratios. 
Several astrophysical mechanisms have been suggested that 
can lead to a number of visible satellite galaxies similar 
to that observed. Perhaps haloes were able to form a few 
stars initially, but the baryonic components of all haloes 
below some critical mass were subsequently destroyed by 
supernova feedback (e.g . iLarsonl 11974 : iDekel fc SUM 11989 : 
iFerrara fc Tolstovl [2000). Alternatively (or perhaps addi- 
tionally) photoionisation may have prevented star formation 
in the smallest haloes (e.g. lEfstathioull 19921 : ISomervillell2002l: 
iHoeft et all l200rj ; ISimon fc Gehal 120071 1. As ISawala et al l 
< |201Ch have shown, these two mechanisms can combine to 
produce very high mass-to- light ratios in haloes of 1O 9 M0 
and below, perhaps reconciling the number of very faint 
dwarf galaxies produced in ACDM simulations with the ob- 
servations. 

In this work, we turn our focus to more massive dwarf 
galaxies, and follow the evolution of the objects that form 
in haloes of 1O 1O M0. Our initial conditions are based on six 
haloes selected from the M illennium-II Simulation (MS-II, 
iBovlan-Kolchin et all 120091 ') . and resimulated at high reso- 
lution using smoothed particle hydrodynamics (SPH). Our 
simulations include cooling and star formation, supernova 
feedback, metal-enrichment and a cosmic UV background. 
Starting at redshift z = 49, we are able to follow the forma- 
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tion of each individual halo and its central galaxy in their 
full cosmological context, all the way to z = 0. 

On the other hand, the large volume of our parent sim- 
ulation allows us to verify that our sample of resimulated 
haloes is representative of haloes of similar mass, and to 
predict a stellar mass - halo mass relation that can be 
tested against observation. With a box size of 137 Mpc 
and a mass resolution of 9.4 x 1O 6 M0, the MS-II has suf- 
ficient dynamic range to capture the statistics of the as- 
sembly of dark matter haloes between 10 9 and 1O 14 M0. By 
comparing its halo/subhalo mass function to t h e observed 
SPSS stellar mass function of lLi fc White] (|2009l ), lGuo et all 
(|2010l ) derived a typical mass-to-light ratio for each halo 
mass. This analysis assumes a monotonic relationship be- 
tween halo mass and galaxy mass with relatively small scat- 
ter, but does not rely on any other assumptions about the 
processes involved in galaxy formation. We use its result 
to test the viability of our simulations and the underlying 
physical model as a description of the formation of "typical" 
ACDM dwarf galaxies. 

The present work constitutes the first direct com- 
parison of high resolution, hydrodynamical simulations of 
individual dwarf galaxies with the observed abundance 
of such objects. We combine the ability to follow star 
formation self-consistently in individual objects with the 
ability to draw conclusions about the general population of 
dwarf galaxies. 

This paper is organised as follows: We begin in Section[2] 
by reviewing the current status of simulations of the forma- 
tion of dwarf galaxies. Section [3] describes the selection of 
haloes for resimulation and the generation of our high res- 
olution initial conditions, while the numerical methods of 
our hydrodynamic simulations are discussed briefly in Sec- 
tion [4] In Section [5j we show results for six haloes of final 
mass 1O 1O M0, and compare the properties of the galaxies to 
previous work, and to observation. In Section|S] we consider 
the predictions of our simulations for the stellar mass - halo 
mass relation and discuss the discrepancy with that inferred 
from comparing the observed stellar mass function to the 
halo abundance in ACDM simulations. We conclude with a 
summary and interpretation of our results in Section [7] 

Unless stated otherwise, where we refer to the mass of a 
galaxy, we mean the stellar mass M* , whereas the mass of a 
halo includes the total dynamical mass enclosed within T200, 
the radius that defines a spherical overdensity 200 times the 
critical density of the universe. When quoting the results for 
our own simulations, we always use physical mass units of 
Mq, assuming h = 0.73. 



2 REVIEW OF PREVIOUS WORK 



Earlier examples of numerical studies of dwarf galaxy for- 
mation and evolution in ~ 10 10 M p) haloes includ e sim- 



ulations by Peluoessv et al. 


l2004h. Stinson et al. ((20071. 


20091). IValcke et al. (l200Sh. 


Mashchenko et al.1 ( 20081) and 


Governato et al.l (hoiOl). 



The first three have investigated the evolution of dwarf 
galaxies embedd ed in dark matter haloes of constant mass. 
iPelupessv et al.l used initial conditions modelled after dwarf 
irregular galaxy DDO 47, set up with a stellar disk of 



Table 1. Results of earlier numerical simulations 



Reference 


M* 


Mtot 


v c or 0* 




[1O 7 M ] 


[1O 9 M ] 


[km s _1 ] 


Pelupessy et al. (2004) 1 


18 


15 


80 


Stinson et al. (2007) 2 


7.86 


5.0 


15.1 


Stinson et al. (2007) 2 


22 


8.6 


20.1 


Stinson et al. (2007) 2 


38.6 


14 


29.9 


Stinson et al. (2009) 3 


1.72 


11 


16.8 


Valcke et al. (2008) 4 


57.9 


1.1 


35.2 


Valcke et al. (2008) 4 


48.8 


1.1 


30.9 


Mashchenko et al. (2008) 5 


1.0 


2.0 




Governato ct al. (2010) 6 


48 


35 


56 


Governato et al. (2010) e 


18 


20 


54 



Notes: Col. 2: Stellar mass, Col. 3: Halo mass (M20o)i Col. 4: 
Maximum rotation velocity 1,6 or 1-D velocity dispersion 2 ' 3,4 . All 
quantities are measured at z = 0, except for Mashchenko et al., 
where the halo mass is at z = 5 and stellar mass at z = 6.2. 
Remarks: 1 Static initial conditions set up to reproduce the dwarf 
irregular galaxy DDO 47; 2,3 Static NFW profiles with initial 
baryon fractions of 10% 2 and 1% 3 ; 4 Runs DH01 and DH02, as- 
suming static, simplified Kuz'min Kutuzov profiles with a = c = 

4 kpc (DH01) and 6 kpc (DH02). 

1.8 x 10 s M and a gas disk of 1.9 x 10 s M Q inside a dark mat- 
ter halo of 1.5 x 10 10 Mq. Thus, the stellar mass to halo mass 
ratio is not a result of their simulation, but was chosen a pri- 
ori. They show that the star formation behaviour in s uch a 
system is consistent with observations. IValcke et al.l stud- 
ied the formation of dwarf elliptical g alaxies assuming cored 
initial dark matter profiles following iDeionghe fc de Zeeuwl 
(|l988l ). and gas at a mass fraction of 17.5%. Cooling, star 
formation and feedback are included in their simulations, 
making their final stellar masses of 4.9 - 5.8 x 1O 8 M for 
haloes of 4.1 x 10 9 Mq a direct prediction of their models. 
IStinson et al.l also assume fixed dark matter profiles in their 
initial conditions, and perform simul ations that inc l ude co ol- 
ing, star formation and feedback. In lStinson et al.l (|2007l ). a 
fixed initial baryon fraction of 10% is assumed, and stellar 
masses of 7.9 x 10 7 - 3.9 x 10 8 M B are pr o duced in haloes of 

5 x 10 9 - 1.4 x 10 10 M Q . In lStinson et al.1 (|2009f ). the baryon 
fraction is varied, and we also include their result with a very 
low initial baryon fractio n of 1%, that lead s to a smaller 
stellar mass IValcke et al.l and IStinson et al.l do not include 
a UV background, which may contribute to the high star 
for mation efficiency in t heir s imula tions. 

iMashchenko et al.l (|2008l ) and iGovernato et al.l (|2010l ) 
both performed simulations that include the formation 
of the dark ma tter halo in a cosmological volume. 
IMashchenko et all used constrained initial conditions, aimed 
at reaching a halo mass of 1O 9 M0 at 2 = 6, and followed the 
evolution up to z = 5. At this time, their halo reached a mass 
of 2 x 1O 9 M0, with 1O 7 M0 of stars formed. By comparison 
with the typical evolution of haloes in the MS-II (see Fig- 
ure^, we note that this is consistent with a mass of 1O 1O M0 
at z = 0. The extent of additional star formation up to z = 
is unknown, however. Most recently, IGovernato et al] have 
performed hydrodynamical simulations of two dwarf irregu- 
lar galaxies at very high resolution, which they follow up to 
2 = 0. This makes these most comparable to our own simula- 
tions, and also makes their results most directly comparable 
to observed, present-day dwarf galaxies. Their simulations 
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Figure 1. Comparison of Halo 4 at z = in a pure dark matter resimulation and in the parent Millcnnium-II Simulation. The left panel 
shows the position of 0.5% of the particles in a box of sidelcngth 1 Mpc in the resimulation, while the central panel shows the position 
of all particles within the same region in the MS-II. The panel on the right shows all particles in a box of 5 Mpc in the MS-II, with 
Halo 4 in the centre. All three panels are centred on the same absolute coordinates for the parent box of sidelcngth 137 Mpc, showing 
the position of the halo to be in perfect agreement. The FoF mass of the halo also agrees to within less than 1%. A comparison of the 
left and the central panel reveals the additional substructure resolved in the resimulation. 



start with values of Sl m = 0.24, Sib = 0.042, and predict 
stellar masses of 1.8 and 4.8 x 10 8 M Q in two haloes of 2.0 
and 3.5 x 1O 1O M0, respectively. 

While these five sets of simulations vary in the setup of 
the initial conditions, the cooling, star formation and feed- 
back recipes, the treatment of the cosmic UV background, 
the simulation code and the numerical resolution, they all 
predict final stellar masses consistent with ~ 10 s Mq for dark 
matter haloes of ~ 1O 1O M0. We give an overview of some of 
the relevant properties of these simulations in Table [T] 



3 INITIAL CONDITIONS 

The parent simulation, as well as our high-resolution resimu- 
lations, are performed in the context of a ACDM cosmology, 
with Q A = 0.75, Sl m = 0.25, h= 0.73 and a s = 0.9, identical 
to the values used f or the original Millennium Simulation 
jSpringel et alj|2005h . 

The Millennium-II Simulation followed structure forma- 
tion in a volume of 137 3 Mpc 3 using 2160 3 dark matter par- 
ticles and periodic boundary conditions. This corresponds 
to a mass resolution of 9.43 x 1O 6 M0, and a force resolution 
of 1.37 kpc. At z = 0, it contains about 12 million friends-of- 
friends (FoF) haloes with at least 20 particles, corresponding 
to a minimum resolved halo mass of 1.9 x 10 8 M o . Haloes 
of ~ 1O 1O M0, the mass at z = that we select for our 
resimulations, are resolved with over 10 3 particles. 

Out of more than 10 4 haloes within our mass range 
in the MS-II, we identified 25 haloes as resimulation candi- 
dates, based on the condition that all particles within twice 
the virial radius at redshift z — were in a connected re- 
gion and inside a sphere of radius 0.67 Mpc in the initial 
conditions. Out of these, six haloes were selected in order 
to study a varied but representative sample of mass accre- 
tion histories (see Figure [2]). The initial conditions for the 
resimulations were generated by re-sampling the region of 
interest with a high number of low mass dark matter parti- 
cles, while the remaining volume was sampled with increas- 
ingly coarser resolution, sufficient to capture the long-range 



tidal field. To account for the higher Nyquist frequency of 
the resimulations, small-scale fluctuations were added to the 
displacement and velocity fields of the original MS-II using 
the method of second-order Lagrangian perturbation theory 
described in IJenkinsI (|2010l ). 

Figure [T] shows the final distribution of dark matter par- 
ticles in slices centred on a dwarf halo in the MS-II and in 
one of our high-resolution resimulations (pure dark matter). 
The difference in mass between the central FoF halo in the 
resimulations and in the MS-II is < 1%, equivalent to ~ 10 
particles in the parent simulation. The position and veloc- 
ity are well reproduced, and the agreement also extends to 
substructures outside the main halo. The structure resolved 
in the MS-II is again found at the correct mass and loca- 
tion, while some additional substructure is visible only in 
the resimulation. 

We have performed resimulations both with pure dark 
matter and with gas particles added, splitting each high- 
resolution dark matter particle at a mass ratio of S7b = 0.046 
to Qom = 0.204. All resimulations start at z = 49 and 
are evolved up to z = 0. The hydrodynamical simulations 
include 1.1 x 10 6 high resolution dark matter particles of 8 x 
1O 4 M0, and an equal number of gas particles of 1.8 x 1O 4 M0. 
The initial stellar particle mass is 9 x 1O 3 M . At z = 0, the 
individual haloes are resolved with more than 10 s particles 
within r2oo- The volume outside of the high-resolution region 
is sampled with an additional 7.6 x 10 5 dark matter particles 
of varying mass to include the evolution of structure on large 
scales. 

To check for possible biases due to our selection method, 
we have compared our candidates for resimulation to the 
total population of similar mass haloes in the MS-II. Fig- 
ure [2] shows as solid lines the merger histories of the six 
selected dark matter haloes in the high resolution simula- 
tions, together with the typical mass accretion history, de- 
rived from merger trees of ~ 10 4 randomly selected, simi- 
lar mass haloes. At each redshift, the inner and outer grey 
regions indicate the 3rd, 16th, 84th and 97th percentiles, 
equivalent to la and 2a deviations from the mean mass for 
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Figure 2. Evolution of FoF-halo mass as a function of redshift 
in our pure dark matter resimulations. The solid coloured lines 
show the mass accretion history of the six haloes we have re- 
simulated at high resolution. Ovcrplotted as a thick dashed line 
is the mean halo mass from the Millennium-II Simulation, for 
all haloes of similar final masses. Also shown are the lcr and 2a 
upper and lower bounds at each redshift, in dark and light shades, 
respectively. It can be seen that there is a variety of assembly 
histories, both in the parent simulation, and in the sample of 
rcsimulatcd haloes. Haloes 3 and 6 have undergone recent major 
mergers, while Halo 5 formed significantly earlier than the five 
others. The variety within our sample is somewhat greater than 
expected for a random sub-sample of the MS-II, but there is no 
systematic bias in formation history. Formation redshifts, defined 
as the time when a halo reaches half its peak mass, lie between 
1 and 2.5, consistent with a median formation redshift of 2. Note 
that the halo mass m2oo can differ from the FoF-mass by ~ 5 — 
15%, and due to outflows, the total mass in each of the six haloes 
is reduced to ~ IO 1O M0 in the resimulations with gas. 



a Gaussian mass distribution. It can be seen that the vari- 
ance within our sample is higher than the expected variance 
within a random sub-sample of haloes. This allows us to fol- 
low the evolution of haloes with a range of merger histories 
in a limited number of simulations. However, there is no sys- 
tematic bias in the mass accretion history of our haloes, so 
our sample can be considered a reasonably unbiased repre- 
sentation of 10 10 Me) haloes in the MS-II. 



galaxies (jSawala et al 



2008 



20091 ) . and of dwarf spheroidal 
2010h . In Sections O to g3J we 
briefly explain the most important characteristics of this 
code, and refer the interested reader to the above references 
for a more detailed description. 



4.1 Gravitational Softening 

In order to reduce two-body interactions arising from the 
particle representation of the matter distribution, the grav- 
itational potential is modified by replacing the divergent 
1/r 2 dependence with l/(r 2 + e 2 ), w here e is the gravita- 
tional softening scale (|Aarsethlll963l ). The choice of e rep- 
resents a compromise between the errors due to residual 
two-body effects, and the loss of spatial resolution below 
several softening scales. We begin our simulations with a 
softening length fixed in comoving coordinates to l/10th of 
the mean interparticle spacing for each particle type, cor- 
responding to ~ lh _1 kpc in the high resolution region. 
After the collapse of the halo, we keep the softening scale 
in this region consta nt in physica l coord inates from z — 7, 
at a value of 155 pc. I Power et all ||2003h give a lower limit 
face = ?~200 / VN200 to prevent strong discreteness effects in 
haloes, which corresponds to ~ 140 pc for a 1O 1O M0 ob- 
ject resolved with A200 ~ 10 5 particles. We also resimulated 
one of our haloes, Halo 4, with a physical softening scale 
of 77.5 pc, and checked that this did not alter the results 
significantly. 



4.2 Cooling and UV Background 

Above the hydrogen ionisation temperature of 10 4 K, our gas 
co oling model is based on metal -dependent cooling functions 
of ISutherland fc Dopital (|l995t ). The model assumes colli- 
sional excitation equilibrium, and does not include metal 
or molecular cooling below 10 4 K. In addition, we include 
Compton cooling, which is the main coolant at high redshift. 
It depends on the free electron density, as well as on the tem- 
perature difference between the gas and the evolving CMB. 
For this purpose, the ionisation states of H, He, and the free 
electron number density are compu ted analytically, follow- 
ing the model of iKatz et ail l|l996ft . We have included UV 
background radiation in our model, which adds a heating 
term to the net cooling function of the partially ionised gas. 
In all simulations, the UV background is present from z — 6, 
and its spectral energy distributi on and the time evolutio n 
of its intensity follow the model of lHaardt fc Madaul (|l996l ). 
A test simulation of Halo 4 without the UV background pro- 
duced over twice as many stars by z = 1, compared to the 
simulation which includes UV radiation. 



4 NUMERICAL METHODS 

The high-resolution simulations presented here h ave been 
perfo r med using the Tree-P M code GADGET-3 (|Springell 
120051 : ISpringel et al.l l2008h . which includes gravity and 
smoothed particle hydrodynamics. As an extension, metal- 
dependent cooling, star formation, chemical enrichment 
and energy injection from type II and type la supernovae 
have been implemented in th e multiphase gas model of 
IScannapieco et all (|2005l , 120061) . This model has previously 



4.3 Star Formation Criteria 

Cold gas particles can spawn, or be converted into, star par- 
ticles, subject to certain conditions. We require the gas par- 
ticle to be in a region of convergent flow. In addition, we im- 
pose a physical density threshold p c on the local gas density. 
The existence of a thr eshold for star fo r matio n is motivated 
by observations (e.g. iKennicuttl 1 19891 . [l998h . Calculations 
bv iQuirkl (Il972h as well as numerical simulations, e.g. by 
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iKatz et all (| 19961 ) : ISpringel fc Hernquistl (|2003h ;[ Bush et al.l 
( 20081 ) and others have shown that the observed Kennicutt- 
Schmidt relation can be reproduced in disk galaxies by 
imposing a volume den sity threshold, even thoug h differ- 
ent values are assumed. iKovama fc Ostriker] (|2009l ) demon- 
strated with high-resolution simulations of the turbulent in- 
terstellar medium that the star formation rate depends only 
weakly on the choice of p c , and valu es in the range 0.1 cm" 3 
(jStinson et all 120091 ) to 100 cm" 3 dGovernato et alJ l2010l ) 
can be found in the recent literature. iGovernato et al.l re- 
ported better convergence in their high-resolution simula- 
tion with a choice of 100 compared to 0.1. In this work, we 
adopt a value of 10 cm -3 . We have also tested a density 
threshold of 0.1 cm -3 , more similar to our own previous 
work. In this case, star formation starts at higher redshift 
and is less bursty. For Halo 4, the final stellar mass increases 
by ~ 36% with a threshold of 0.1 cm" 3 . This difference is less 
than the variance in stellar mass between individual haloes, 
and does not qualitatively affect the stellar mass-halo mass 
ratio. The limited effect of p c results from the fact that star 
formation is mostly self-regulating in our simulations. We 
also impose a threshold p g /pg 10 4 on the local gas over- 
density relative to the cosmic mean, which ensures that star 
formation only takes place in virialized regions even at very 
high redshift. 



4.4 Star Formation Efficiency 

Subject to the constraints described in Section [4.31 the star 
formation efficiency is regulated by a single efficiency param- 
eter c*, so that the star formation rate density is given by 
p* = c *Pgt^y n , where td yn is the local gas dynamical time. 
The creation of an individual stellar particle of mass m* 
from a gas particle of mass m g during the time interval At 
is stochastic, with probability given by: 



P* = 



exp 



At 

Idyn 



In simulations with radiative transfer, iGnedin et al.l 
(|2009l ) found that dust acts as a catalyst for molecular cloud 
formation, suggesting that star formation may be less ef- 
ficient in the low metallicity environment of dwarf galax- 
ies. As our simulations cannot follow cloud formation, we 
do not take this in to account and assume a constant c*. 
iRicotti et al. I l|2002T l showed that if star formation is strongly 
self-regulating, the star formation rate is determined primar- 
ily by the thermodynamic properties of the gas and depends 
only very weak ly on c*. This resu lt was confirmed in our 
previous work (Sawala et al.ll20ld ). and in all simulations 
presented here, we adopt our earlier value of c* = 0.05. 

Each star particle is produced with a single stellar 
population, whose metallicity is inherited from the parent 
gas particle. W e assume a Salpeter initial mass function 
|SalpetedlT955T ), and calculate stellar luminosities using the 
stellar synthesis model of iBruzual fc Charlotl (|2003l ). 



4.5 Multiphase Interstellar Medium and Feedback 

For each star particle, we determine the rate as well as 
the yields of supernovae type II and type la. Chemical 



yie lds are calculated separat ely for the two types, follow - 
ing lWooslev fc Weaver! (|l995l ) and I Thielemann et alJ l|l993l ). 
respectively. Supernovae type II are assumed to be instan- 
taneous, while supernovae type la follow a uniform delay 
time distribution between 100 Myrs and 1 Gyr. We assume 
a constant energy production of 7 x 10 50 ergs per super- 
nova, which is released into the interstellar medium (ISM) 

as thermal energy. 

The multiphase scheme of lScannapieco et all (|2006l ) al- 
lows an overlap of diffuse and dense gaseous components. 
This preserves the multiphase structure characteristic of the 
ISM, in which components with very different temperatures 
and densities coexist. It also avoids the overestimation of 
density in diffuse gas near high density regions which can 
cause a serious underestimate of its cooling time. The de- 
coupling is achieved by considering as neighbours in the SPH 
smoothing kernel only gas particles with similar thermody- 
namical properties, as defined by the ratio of their entropic 
functions. 

Supernova energy is shared equally between the hot and 
cold phases. Cold particles which receive supernova feedback 
accumulate energy until their thermodynamic properties are 
comparable to those of their local hot neighbours. At this 
point, the energy is released and t he particles are promoted 
to the hot phase. We supplement IScannapieco et all (|2006l ) 
by including a seeding mechanism that defines reasonable 
properties for the local hot phase, even if no neighbouring 
particles are considered hot at the time. A cold gas parti- 
cle which has received sufficient supernova energy to raise 
its thermodynamic properties to this level can thus be pro- 
moted, even if it currently has no hot neighbours. This en- 
sures that the distribution of supernova feedback is not de- 
layed at the earliest stages of star formation, when the entire 
interstellar medium can be in a cold and dense configuration 
(see Figure [3]). The seeding mechanism does not create heat 
artificially and conserves energy. We have checked that the 
amount of gas required to seed the hot phase is small, and 
that the ensuing evolution of the two phases is consistent. 



5 GALAXY FORMATION 

The different merger histories of the dark matter haloes de- 
scribed in Section [3] are reflected in their gas accretion his- 
tories, and in the evolution of the galaxies that form within 
them. Section \5. 1 1 describes the co-evolution of the halo and 
its galaxy, while Section 15.21 discusses the properties of the 
final objects. 



5.1 Galaxy Evolution 

The halo assembly histories vary significantly, both among 
our selected sample and among the total population of 
1O 1O M0 haloes. As shown in Figure [21 haloes 3 and 6 had 
recent major mergers at z — 0.32 and z = 0.21, respectively, 
whereas haloes 1 and 2 experienced their last major mergers 
around z = 1. Haloes 4 and 5 have not undergone any major 
mergers since before z = 2. Halo 5 is also significantly more 
massive compared to the other haloes at high redshift. 

The six panels in Figure show the gas mass bound 
to each of the six haloes as a function of redshift. While 
the total coloured area indicates the total amount of gas 
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Figure 3. Evolution of total gas mass and star formation rate as a function of redshift. The top row shows, from left to right, haloes 1-3, 
while the bottom row shows haloes 4-6. The coloured area indicates the total amount of gas in each halo in units of Mq, corresponding 
to the scale on the left. The colour coding indicates the differential amount of gas at a given mass density in units of g cm -3 , as denoted 
by the colour bars above. Overplotted in black is the star formation rate as a function of redshift, in units of Mq yr" 1 , corresponding to 
the scale on the right. In each halo, star formation is tightly coupled to the amount of dense gas, and occurs in bursts, often separated 
by several 100 Myrs, and associated both with supernova feedback and with mergers. Halo 5 assembles earlier than the five others, 
and its significantly higher mass at early times leads to a more prolonged starburst in the galaxy contained within it. All galaxies are 
star-forming at z = 0. 



in each halo, the star formation rate in the main pro- 
genitor, overplotted in black, depends on the presence of 
cold and dense gas, shown in red and orange colours. We 
find that two different mechanisms lead to a burstiness of 
star formation in our simulations. On timescales of hun- 
dreds of Myrs, self-regulation of star formation and super- 
nova feedback lead to periodic variations in the gas den- 
sity, and periodic st ar formation behaviour . This confirms 
the ea r lier results of Pelupessv et al.l ||2004). Stinson et al.l 
(120071). [Mashchenko et all <|2008h . IValcke et all (|200Sl ) and 
iRevaz et all (|2009h . In addition, gas-rich mergers can induce 
starbursts. These bursts are irregular and can be separated 
by several Gyrs. For example, the star formation episodes in 
Halo 4 beginning at z = 8, 4.3 and 2.1 are preceded by merg- 
ers at z = 8.5, 4.7 and 2.4 with haloes of 2 x 10 s , 3 x 10 s and 
10 9 Mq, respectively, which bring in fresh gas. Observations 
of periodic bursts lasting several hundred m illion years have 
been reported for three dwarf galaxies by iMcQuirm et al 
1120091 ). while a number of dw arf gala xies (Leo I (jDolph: 



20021), Leo A (ICole et alj|2007t ). I C 10 dCole!l201Ch. IC 1613 
( Skillman et al.ll2003l). DDO 210 (|McConnachie et al.ll2006l ) 
and Carina (|Koch et al.ll2006l )) show ex tended, quie scent pe- 
riods between star formation epochs. IColel (|201Cf ) suggest 
mergers and gas accretion as triggers for star formation, but 
note that individual bursts and mergers can no longer be 
linked observationally after several Gyrs. 

The star formation history of each individual galaxy 



in our simulations reflects a combination of internal self- 
regulation via supernova feedback, and the supply of fresh 
gas via accretion and mergers. These two effects largely de- 
termine the variance in stellar mass between the haloes in 
our simulations; while differences in merger histories increase 
the variance, self-regulation via feedback decreases it. In our 
sample of six haloes of equal final mass, the galaxy stellar 
masses vary by about a factor of two. 

Gas-rich mergers after z = 6 imply that the progeni- 
tors did not lose all their gas due to the UV background. In 
our simulations, such mergers occur wi th haloes that would 
reach masses above ~ 10 9 M Q by z = 0. ISawala et~aT1 (|20ld ) 
showed that at this mass, a combination of UV and super- 
nova feedback removes gas efficiently, while UV radiation 
alone is not always s ufficient. Observatio ns of local group 
dwarf spheroidals (e.g. iMonelli et al.|[201ol ) also suggest that 
reionisation had at most a minor effect on these galaxies. 

Some major mergers also contribute stars. The fraction 
of final stellar mass formed outside of the main progenitor 
ranges from ~ 5% in Halo 5, accreted at z = 5.2, to close to 
40% for Halo 1, resulting from two major mergers at z — 3.4 
and z = 0.8. Haloes 2 and 6 both accrete ~ 25% in mergers 
at z = 0.4 and 1, respectively, while haloes 3 and 4 accrete 
~ 7% in mergers at z = 0.5 and z = 2.4. Since haloes 3, 4 
and 5 follow a more typical assembly history, we expect the 
typical fraction of stars formed outside the main progenitor 
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Table 2. Overview of numerical simulation results 



Halo 


M* 


M g 


M DM 


ri/2 


0* 




Zmm 




[1O 7 M ] 


[10 7 M o ] 


[1O°M ] 


[kpc] 


[km s _1 ] 


[km s — "'"kpc] 




1 


7.81 


2.81 


9.41 


0.87 


21 


1.4 


0.77 


2 


5.25 


6.23 


8.34 


0.39 


17 


0.7 


0.96 


3 


4.94 


19.4 


8.80 


0.28 


15 


1.1 


2.65 


4 


8.14 


4.17 


8.98 


0.36 


19 


0.1 


2.37 


5 


10.2 


7.11 


9.94 


0.74 


33 


19 


5.17 


6 


6.17 


5.07 


8.54 


0.68 


18 


0.3 


0.38 



Notes: Col. 2: Stellar mass, Col. 3: Gas mass, Col. 4: Dark matter 
mass (all measured inside T200 at z = 0), Col. 5: Stellar half-mass 
radius, Col. 6: De-projected 1-D RMS stellar velocity dispersion, 
Col. 7: Specific stellar angular momentum L = |r X v (<r* and 
L* measured within 3 kpc), Col. 8: Rcdshift of last major merger 
(progenitor mass ratio < 3:1). 

in dwarf galaxies of M* ~ 1O 8 M0 to be < 10%, albeit with 
possible exceptions. 

In general, we find that the mean metallicity evolves 
with age, indicating the recycling of enriched gas in subse- 
quent generations of stars. At each stellar age, we also find 
a spread in metallicities, which indicates the incorporation 
of fresh material. However, we note that due to a lack of dif- 
fusive metal mixing, the metallicity spread in our simulated 
stellar populations can be as high as 3 dex, which is larger 
than observed. 



5.2 Galaxy Properties 

Several properties of the six simulations are listed in Table[2] 
They appear to be in broad agreement with the previous 
studies discussed in Section [2] The final baryon fraction of 
the haloes are between 1.1 and 2.7 % of the total matter. The 
final stellar masses of the six galaxies fall between 4.9 x 10 7 
and 1.0 x 1O 8 M0, which corresponds to stellar mass to total 
mass ratios in the range of ~ 5 x 10" 3 - 10" 2 . All galaxies 
are star-forming at z — 0. 

The final gas masses vary from 2.8 x 10 7 to 1.9 x 10 s Mq. 
Observed dwarf galaxies of this stellar mass typically have 
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a substantial HI content (e.g. 

iGeha et al]|200rj ; Ide Blok et al.l I2008I V In our simulations, 
we cannot directly measure the amount of HI gas. Defin- 
ing the cold gas as the total amount of gas at temperatures 
below the peak of the cooling curve and correcting for the 
contribution of helium and metals, we derive approximate 
upper limits for the HI m asses between lOjand 1.2 x 10 s M . 
We adopt the notation of IGeha et alj (|2006) in defining the 
HI mass fraction as /hi = Mm /(Mm + M,), noting that 
He, H2, hot gas and metals are neglected in the denomina- 
tor. We derive upper limits for /hi between 11% and 71% in 
our six galaxies with a median of 35%. In a sample of 101 
flux-selected SDSS dwarf galaxies of similar stellar mass, 
IGeha et al.l find a higher mean HI fraction of 60%. However, 
individual galaxies show a large scatter in /hi, with several 
as high as 95%, and others with upper limits below 10%. 

The range of stellar half-mass radii in our six simula- 
tions is 0.28 to 0.87 kpc. This is comp arable to the typi- 
cal values obtained bv lGeha et al.l (|2006t ) for half-light radii 
(r-band), even though the observed sample also contains a 
handful with half-light radii greater than 1.5 kpc. We note 



that only Halo 5, which has the most quiescent assembly 
history, contains a galaxy with a rotationally supported stel- 
lar disk, as reflected by the specific angular momentum L*. 
The other five haloes have more ellipsoidal morphologies 
an d very little rotation. This is in contrast to the result 
of lGovernato et"aTI |2010I ). who report pristine disk galax- 
ies in both of their simulations, albeit in haloes of 2.0 and 
3.5 x 1O 1O M0. Observations also suggest that isolated dwarf 
gala xies of this stellar mass are more frequently disk-like 
(e.g. iHunter fc Elmegreenll2006l : IGeha et al.ll2006f ). The lat- 
ter find that 30% of edge-on dwarf galaxies show coherent 
rotation profiles, a number which drops to 18% when all axis- 
ratios are included, but this should b e considered as a lower 
limit. In a volume-limited sample, ISanchez-Janssen et al.l 
(|2010l ) report that less than 30% of galaxies with stellar 
masses of 1O 8 M0 have apparent axis ratios below 0.5, and 
less than 5% below 0.3. They attribute this flattening to 
stellar feedback. Compared to these observations, our simu- 
lated sample of six haloes is too small to assess the statistical 
significance of a 1 in 6 result. It should also be noted that 
the stellar half-mass radii of our galaxies are only resolved 
with a few softening lengths. Therefore, results on details of 
the galactic structure from our simulations are inconclusive. 



6 STELLAR MASS - HALO MASS RELATION 

In Section [3l we explained how we constructed our initial 
conditions from the parent simulation. Its large volume and 
high dynamic range allows us to demonstrate that we have 
resimulated a representative sample of haloes. We can there- 
fore derive implications for the global population of galaxies 
that form in similar mass haloes, and compare with expec- 
tations from matching the observed stellar mass function to 
the abundance of haloes in a ACDM universe. 

Assuming a monotonic relat ionship be t ween stellar 
mass and maximum halo mass, iGuo et all (|2010T ) have 
compared the abundance of haloes/subhaloes in the Mil- 
lennium and Millennium-II Simulations to the observed 
abundance of galaxies as a function of stellar mass obtained 
from the SDSS DR-7 by iLi fc White] (|2009l ). The observa- 
tional sample contains over half a million galaxies at low 
redshift, and extends down to stellar masses of 1O 83 M 
with very small error bars. The combination of the two 
very large simulations also leads to very small errors on 
the theoretical halo abundance. The derived stellar mass 
to halo mass ratio peaks at Mhaio = 1O 11,8 M0 at a star 
formation efficiency of about 20%, and decreases both 
for more massive and for less massive haloes (see Figure 
[4|. The decrease in efficiency at the high mass end is 
generally attributed to AGN feedback (e.g. ICroton et al.l 
I2006J : iBower et al.l 120071 ) . while the decrease for lower mass 
haloes is assumed to be due to the increasing efficiency of 
supernova feedback, and the effect of the UV background. 
This general be haviour was noted e arlier from lower 
recis io n data by Navarro fc Steinmetz J 20001). Yang et "all 
20031), iDekel fc Wool (120031 ), IConrov fc Wechslerl (|2009T l 
and lMoster et all |2010r i. 



Following lYang et all (|2003f ), lGuo et alj (|201Ct) adopted 
the following functional form for the mean stellar mass to 
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They report an accurate fit to the data with parameters 
c = 0.129, Mo = 10 114 M Q , a = 0.926, /3 = 0.261 and 
7 = 2.440. At the low mass end, the SDSS DR-7 data of 
iLi fc White! extends to stellar masses of 2 x 10 8 Mq with 
high accur acy, corresponding to a halo mass of 1O 1O ' 8 M . 
IGuo et al] extrapolate this relation down to the lowest halo 
masses resolved in the MS-II. For halo masses of ~ 10 10 M Q , 
this predicts a stellar mass of ~ 8 x 1O 5 M . 

From a similar analysis based on S DSS DR-3, wh i ch ex - 
tends to stellar masses of 3.2 x 10 8 Mm. IMoster et al.l (|201Ch 
also derive a stellar mass - h alo mass relation, in good agree- 
ment with lGuo et al.l (|2010t ) at the high mass end. They pre- 
dict haloes of 10 10 Mq to host galaxies with stell ar masses of 
~ 5.7xl0 6 M Q , significantly higher than found bv lGuo et al.l . 
but still an order of magnit ude lower tha n our hydrody- 
namical simulations predict. iMoster et al.l note, however, 
that such haloes are at best marginally resolved at their 
mass resolution (2.8 x 10 8 Mq), prohibiting a self-consistent 
treatment of subhaloes. Consequently, they only include 
haloes above 1.6 x 1O 1O M0 in the conditional mass function. 
IMoster et al.l also apply their analy sis to an analytic Sheth- 
Tormen mass function obtained by IVale fc Ostrikerl (|2006h . 
In this non-parametric model, a halo mass of 1O 1O M0 cor- 
responds to a stellar m ass of 1.9 x 1O 6 M , more similar to 
the value of lGuo et all (|201Ch . 

In Figure [ 4] we p lot the stellar mass - halo mass relation 



of IGuo et all (|2010h for haloes between 10 9 and 1O 12 M . 
The solid section of the line shows the relation in the re- 
gion directly derived from SDSS DR-7 data where the un- 
certainties are very small. The dashed section denotes an 
extrapolation to stellar masses below 1O 83 M , assuming a 
faint-end slope of a = —1.15 for the stellar mass function, 
as reported by ILi fc White! HH). Studies of the faint-end 
of the stellar mass function are either limited to nearby re- 
gions or galaxy clusters, or require corrections for incom- 
pleteness and, in the case of photometric redshifts, back- 
ground s ubtraction, which intr oduce considerable uncertain- 
ties fe.g. IChristlein et al.ll2009l ). As a result, different values 
for a in the ra nge of —1.1 to —1.6 are found in the recent 
literature (e.g. iTrentham et al.l |2005l: jBlanton et all 120051 ; 
ICarrasco et al.1120061 : IBaldrv et al.ll2008h The dark grey area 
in Figure [4] shows the effect of a steepening of the faint-end 
slope up to a = —1.58, the value reported bv lBaldrv et al.l 
(120080 . While this has a strong effect on the lowest mass 
haloes, we note that it cannot account for the discrepancy 
we find in haloes of 1O 1O M . In order to fit the constraints 
of SDSS DR-7, the maximal dispersion at fixed halo mass is 
0.2 dex in M*, indicated by the light-grey area. We overplot 
the results of our six simulations as red squares and add 
other 2 = predictions from the studies listed in Table [JJ 
correcting all halo masses for baryonic effects, as described 
below. It is apparent that all these hydrodynamical simula- 
tions overproduce stellar mass for their respective halo mass 
by at least an order of magnitude. 

In Table [3] we compare the properties of our six simula- 
tions to the abundance matching predictions. We note that, 
due to the outflow of baryons, the total mass of our six haloes 
is almost a factor of 1 — fii,/fi m smaller than the masses of 
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Figure 4. The stellar mass - halo mass relation, derived by Guo 
et al. (2010), compared to the results of several numerical sim- 
ulations. The solid black line denotes the range constrained by 
the SDSS DR-7 data, while the dashed line is an extrapolation to 
lower masses, with a faint-end slope of —1.15 for the stellar mass 
function. The dark grey area shows the influence of a steeper 
faint-end slope up to —1.58 (dash-dotted line), while the light 
grey area represents the maximally allowed dispersion of 0.2 dex 
in M* for a given halo mass. The coloured symbols denote the 
results of hydrodynamical simulations, as listed in Table [JJ and 
Table |21 excluding simulations that did not evolve to z = or as- 
sume a baryon fraction of 1% ab initio. The red squares indicate 
our own six simulations. For consistency, we use the halo masses 
of the pure dark matter simulations for our own simulations, and 
apply a correction of 20% to all other haloes, where this informa- 
tion is not available. All hydrodynamical simulations overprcdict 
the stellar mass with respect to the observed relation by more 
than an order of magnitude. 



the corresponding haloes from the pure dark matter sim- 
ulation. This effect is expected at such low star formation 
efficiency. For consistency with lGuo et al.l (|2010l ). we there- 
fore use the (higher) peak masses of the pure dark matter 
simulations in deriving the stellar mass predicted for each 
of our haloes by the abundance matching argument. For all 
galaxies listed in Table [JJ whose peak halo mass cannot be 
defined or is not given, we increase the halo mass in Figure|4] 
by Qt/Qm ~ 20%, the maximally expected correction. 

Comparing the results of our simulations to the predic- 
tions, we find that the hydrodynamical simulations overpro- 
duce stellar mass by a median factor of ~ 50. Alternatively, 
abundance matching predicts that galaxies of 1O 7 ' 9 M , the 
median stellar mass produced in our hydrodynamical sim- 
ulations, should reside in haloes with typical masses of 
~ 4.5 x 1O 1O M , rather than 1O 1O M . If 1O 1O M haloes re- 
ally hosted galaxies with M* = 1O 7 9 M , a ACDM universe 
would overpredict their abundance by a factor of ~ 4. 

This discrepancy is too large to be attributed solely 
to incomplet e ness in the observed stellar mass function. 
iBaldrv et al. (2008) have used the stellar mass - surface 
brightness relation of SDSS galaxies in order to estimate 
the complet e ness at the faint end. Based on this analysis, 
ILi fc White] (|2009l ) estimate the completeness at 1O 83 M 
to be well above 70%. Following IBaldrv et all the uncer- 
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Table 3 


Comparison of stellar 


mass — 


halo mass ratios 


Halo 


M* 


Mtot 


M max 


M^(SMF) 




[1O 7 M ] 


[1O 9 M ] [ 


1O 9 M ] 


[1O 7 M ] 


1 


7.81 


9.52 


12.1 


0.154 


2 


5.25 


8.46 


11.2 


0.1216 


3 


4.94 


9.04 


10.7 


0.104 


1 


8.14 


9.10 


11.8 


0.143 


5 


10.2 


10.1 


12.5 


0.171 


6 


6.16 


8.65 


10.7 


0.104 



Notes: Col. 2: Stellar mass obtained from simulation, Col. 3: Com- 
bined mass (111200) °f stars, gas & dark matter of the halo in the 
hydrodynamical simulation, Col. 4: Peak halo mass (ni2oo) m the 
pure dark matter simulation, Col. 5: Stellar mass corresponding 
to M max from the abundance matching. 

tainty in the number of 1O 7 ' 9 M0 galaxies is much smaller 
than the discrepancy we report. 

The difference is also unlikely to be attributable to nu- 
merical errors in our hydrodynamical simulations, or to the 
specific parametrisation of star formation and feedback in 
our model. From Table[TJ it is clear that all other current hy- 
drodynamical models, while succeeding in reproducing many 
of the observed features of individual dwarf galaxies, also 
predict similar or higher galaxy mass to halo mass ratios. 
They thus also fail to reproduce the low star formation ef- 
ficiencies required to explain the observed abundances of 
dwarf galaxies in a ACDM universe. 

Haloes of 10 10 M o are well resolved in the Millennium- 
II Simulation, and the number of such haloes in a volume 
of 137 3 Mpc 3 is clearly large enough for statistical uncer- 
tainties to be small. While the assumed cosmological pa- 
rameters of = 0.75, fi m = 0.25 and <j% — 0.9 are 
only marginally consi stent with the five-year WMAP data 
|Komatsu et alj 120091 ). this has a negligible effect on the 
mass function. It is unlikely that the number density of 
1O 1O M0 haloes formed in a ACDM cosmology is significantly 
overestimated in our parent simulation. A lower value of erg 
would slightly inc r ease t he abundance of haloes of this mass, 
but as I Yang et all (|2003h have shown, the abundance match- 
ing result, which depends on the cumulative abundance of 
more massive haloes, is almost unchanged at 1O 1O M0 in 
ACDM. 

In warm dark matter (WDM) models, structure is 
erased below a characteristic free-streaming length that 
depends on t he a ssumed properties of WDM particles. 
IZavala et al.l (|2009T ) have compared the halo mass func- 
tions in high resolution dark matter simulations of ACDM 
and AWDM, assuming ttiwdm = 1 keV, and find that the 
present-day abundance of haloes of 1O 1O M0 decreases by a 
factor of ~ 3. These simulations truncate the power spec- 
trum of the initial conditions in the WDM case but neglect 
thermal velocities, which could increase the effect. How- 
ever, recent combined analysis of structur es observed in the 
Lyman-Q forest of SDSS and HIRES by IViel et all (^OOct ) 
suggest a lower limit of ~ 4 keV for thermal relics in a pure 
AWDM cosmology, which would allow only a much smaller 
deviation from CDM on these scales. While a WDM model 
could thus perhaps account for the reported discrepancy 
between simulations and star formation efficiencies inferred 
from abundance matching, the required WDM particle mass 



appears disfavoured by observation. A WDM solution would 
also significantly alter the internal structure of dwarf haloes, 
and simultaneously reduce the abundance of lower mass ob- 
jects. 

In principle, one could turn to direct measurements of 
halo masses for individual galaxies on the relevant scales, 
to elucidate whether the halo mass function of the ACDM 
model, and the inferred stellar mass - halo mass relation- 
ship, are correct. Direct mass estimates, t hrough gravita- 
tional lensing (e.g. iMandelbaum et~aL I I2006J ) are only avail- 
able for haloes with masses above ~ 1O 118 M0, however, 
where the y agree wit h the CDM predictions and the rela- 
tionship of lGuo et all F or dwarf galaxies, on e has to rely on 
HI rotation curves (e.g. Ide Blok et al.l|2008T ). which do not 
give reliable estimates for total halo masses. 

As demonstrated in Section [3] we have been careful to 
resimulate a representative sample of haloes, and to exclude 
any systematic bias. Considering the limited variance in stel- 
lar mass among the six haloes, statistical fluctuations are an 
unlikely source for the discrepancy. 

It is worth noting that semi-analytical models of ga laxy 
formation (e.g. iKauffmann et al1ll993l ; ICole et al]|2000l ) at- 
tempt to reproduce the observed faint-end slope of the stellar 
mass function in a ACDM universe by assumi ng highly effi- 
cient supernova feedback in small haloes (e.g. iBenson et al.l 
120031 ; iKhochfar et~all 120071). We h ave applied the semi- 
ana lytic model of Guo et all (|2010l . b), which reproduces 
the lLi fc White] (|2009T l stellar mass function, to the merger 
trees of our six resimulated haloes, and also to a randomly 
selected sample of similar mass haloes from the Millennium- 
II Simulation. We find no difference in the predicted stellar 
mass between the selected haloes and the random sample, 
but stellar masses that are roughly two orders of magnitude 
smaller than in our hydrodynamical simulations. Indepen- 
dent of the cause of the discrepancy between the hydrody- 
namical simulations and the observed stellar mass function, 
there is thus a divergence between current hydrodynamical 
and semi-analytical models for dwarf galaxies. The semi- 
analytical models are tuned to produce the correct galaxy 
abundance, whereas hydrodynamical models aim at repro- 
ducing the physical processes and the structure of individual 
galaxies. Clearly, these two aspects cannot be treated sepa- 
rately, if we are to converge to a consistent picture of galaxy 
formation. 



7 SUMMARY 

We have performed high-resolution hydrodynamical simula- 
tions of six ~ 1O 1O M0 haloes, extracted from a large, cosmo- 
logical parent simulation. We find that differences in merger 
histories lead to the formation of dwarf galaxies with differ- 
ent star formation histories and final stellar masses between 
4.9 x 10 7 and 1O 8 M . These stellar masses agree with previ- 
ous simulations of similar mass haloes, and the structure of 
our simulated galaxies resembles that of observed galaxies 
of similar stellar mass, to the extent which we can resolve 
structure in our simulations. 

However, all these simulations imply an efficiency of 
conversion of baryons into stars which is at least an order 
of magnitude larger than that which is required to explain 
the observed abundance of dwarf galaxies in a ACDM uni- 
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verse. While current hydrodynamical simulations, including 
our own, are consistent with almost arbitrarily high mass- 
to-light ratios for the faintest galaxies in haloes of 10 9 Mq or 
less, they thus appear to be inconsistent with the mass-to- 
light ratios of larger dwarf galaxies, even when a moderately 
steep faint-end slope of the stellar mass function is assumed. 
The current recipes for mechanisms such as UV heating and 
supernova feedback appear sufficient to remove the "Miss- 
ing Satellites Problem" for the smallest satellites. However, 
isolated dwarf galaxies with stellar masses of 10 8 M Q are still 
substantially overproduced in current hydrodynamical sim- 
ulations, even when these mechanisms are included. 

Our results suggest three possible explanations: The 
current observational count of dwarf galaxies could be 
incomplete, underestimating the true number density of 
10 8 Mq galaxies by a factor of four or more. In that case, the 
hydrodynamical simulations could be correct, but the semi- 
analytical models that produce low abundances of dwarf 
galaxies have been tuned to incorrect data. 

If the count of dwarf galaxies is almost complete at 
10 8 Mq, these galaxies must, in a ACDM universe, be re- 
siding in haloes significantly more massive than 10 10 M Q , 
and all current hydrodynamical simulations overpredict the 
efficiency of star formation by more than a factor of ten. 
This could be an indication of numerical problems, or, more 
likely, of incorrect or incomplete assumptions about the rel- 
evant astrophysics. Several possible mechanisms may con- 
tribute to a star formation efficiency in current simulations 
that is too high compared to real galaxies: 

• Supernova feedback may be more efficient in ejecting 
gas from dwarf galaxies than current hydrodyn amical sim- 
ulations predict. For example, iGuo et all (|2010l . b) showed 
that the observed stellar mass function can be reproduced in 
semi-analytic models by assuming very strong mass-loading 
of winds in low mass haloes. 

• The full effect of reionisation on the IGM may not 
be captured in current models. As a result, cooling times 
may be underestimated, and the fraction of gas-rich merg- 
ers overestimated. Local sources of extreme UV and soft 
X-ray radiation may also ionise the interstellar medium, in- 
ducing another self-regulation mecha nism for star formation 
i|Ricotti et alj|2002l : ICantalupoll2010l '). 

• Low dust content may lead to less efficient cloud forma- 
tion and shielding at low metallicities. The transformation 
rate of cold gas into stars, currently assumed to be universal, 
may therefore be ov erestimated in dwarf galaxy simulations 
(|Gnedin et al.ll2009T l. 

• Processes such as magnetic fields, cosmic rays and the 
feedback from population-Ill stars are not included in any 
of the current models, and may further reduce the star for- 
mation efficiency. 

Any revised model would however still have to reproduce 
features of individual galaxies consistent with observations. 
We also note that a model which substantially decreases the 
number of 10 s M© galaxies would imply that the halo masses 
of fainter dwarf galaxies would need to be revised upwards, 
as some of these would now be required to live in 10 10 Mo 
haloes. 

If the observed stellar mass function is complete, and 
the hydrodynamical simulations correctly capture the rele- 
vant physics of galaxy formation, the Millennium-II Sim- 



ulation (and similar ACDM simulations) overpredict the 
number of 10 10 Mq dark matter haloes. This would seem to 
require the underlying physical assumptions of the ACDM 
model to be revised. Warm Dark Matter may offer a pos- 
sibility, but only for particle masses of ~ 1 keV, below the 
limit apparently implied by recent Lyman-a observations. 

Of the three proposed scenarios, it appears that missing 
astrophysical effects in the simulations are the most likely 
cause of the discrepancy, and the most promising target in 
search of its resolution. While the three scenarios differ in 
nature, none is without significant implications for galaxy 
formation, which will have to be addressed in the future. 
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